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J.—S ^ We consider wireless mesh networks and the problem of routing end-to- 

^sj ■ end traffic over multiple paths for the same origin-destination pair with 

minimal interference. We introduce a heuristic for path determination 
with two distinguishing characteristics. First, it works by refining an ex- 
tant set of paths, determined previously by a single- or multi-path routing 
K^ ' algorithm. Second, it is totally local, in the sense that it can be run by 

5_j , each of the origins on information that is available no farther than the 

Cu ' node's immediate neighborhood. We have conducted extensive computa- 

tional experiments with the new heuristic, using AODV and OLSR, as well 
as their multi-path variants, as underlying routing methods. For two dif- 
ferent CSMA settings (as implemented by 802.11) and one TDMA setting 
running a path-oriented link scheduling algorithm, we have demonstrated 
that the new heuristic is capable of improving the average throughput 
network-wide. When working from the paths generated by the multi-path 
routing algorithms, the heuristic is also capable to provide a more evenly 
distributed traffic pattern. 
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1 Introduction 

Wireless mesh networks (WMNs) have lately been recognized as having great po- 
tential to provide the necessary networking infrastructure for communities and 
companies, as well as to help address the problem of providing last-mile connec- 
tions to the Internet [331133]. However, mutual radio interference among the net- 
work's nodes can easily reduce the throughput as network density grows above 
a certain threshold iTj and therefore compromise the entire endeavor. Such in- 
terference is caused by the attempted concomitant communication among nodes 
of the same network and constitutes the most common cause of the network's 
throughput's falling short of being satisfactory (hardly reaching a fraction of 
that of a wired network [21] )• A promising approach to tackle the reduction of 
mutual interference seems to be to combine routing algorithms with some inter- 
ference avoidance approach, such as power control, link scheduling, or the use 
of multi-channel radios [2]. In fact, this type of network interference problem 
has been addressed by a considerable number of different strategies to be found 
in the literature [12l[T3l[Il|4ll[IIl[5ll|44l|4]. 

An alternative approach that presents itself naturally is the use of multi- 
path routing to distribute traffic among multiple paths sharing the same origin 
and the same destination, since in principle it can help with both path recov- 
ery and load balancing better than the use of single-path strategies. It may, 
in addition, lead to better throughput values over the entire network [26j [5] . 
But while these benefits accrue only insofar as they relate to how the multi- 
ple paths interfere with one another [JHl Ei] ; unfortunately this aspect of the 
problem is not commonly addressed by multi-path strategies. What happens as 
a consequence is that, though promising by virtue of adopting multiple paths 
to accommodate the same end-to-end traffic, in general such strategies fail to 
perform as desired because they do not tackle the interference problem during 
path discovery. The single noteworthy exception here seems to be the algorithm 
reported in [49) , but it uses geographic information (like localization aided by 
GPS) to find paths with sufficient spatial separation so as not to interfere with 
one another. In our view this weakens the approach somewhat, since such type 
of information may not always be available |14j . Moreover, the corresponding 
algorithm relies on the solution of an NP-hard problem on an input that has 
the size of the network [50], so the solution may be unattainable in practice. 

Here we propose a different approach to alleviate the effects of interference 
in multi-path routing. Our approach is based on two general principles. First, 
that it is to work as a refinement phase over existing routing algorithms, thereby 
inherently preserving, to the fullest possible extent, the advantages of any given 
routing method. Second, that it is to rely only on information that is locally 
available to the common origin of any given set of multiple paths leading to 
the same destination. That is, only information that the origin can obtain by 
communicating with its direct neighbors in the WMN should be used. One 
intended consequence of the latter, in particular, is that refining the set of 
paths departing from any common origin should be easily implementable by 
straightforward message passing, and moreover, that any required calculation 



by that node should be amenable to being carried out efficiently even if it 
involves the solution of a computationally difficult problem. In order to comply 
with these two principles, our approach operates on a previously established set 
of paths leading from a common origin, say i, to a common destination, say j. 
It operates exclusively on the neighborhood information stored at node i itself 
or at any of its neighbors, say k, such that k participates in some of the i-toj 
paths, as well as on the information stored at these same nodes regarding the 
routing of packets to node j. Once node i has acquired all this information, an 
undirected graph G^ is constructed that represents every possible interference 
that can occur as packets get forwarded toward j by those of i's neighbors that 
are on i-to-j paths. Solving a well-known NP-hard problem (that of finding a 
maximum weighted independent set) on this typically small graph serves as a 
heuristic to decide which of the i-to-j paths to keep and which to discard. 

It is important to note that, being determined with reference to graph Gy, 
the resulting set contains no two paths that interfere with each other as far as 
node i's neighbors are concerned, except of course for the inevitable interference 
that may occur as packets leave i or reach j. In our view, this provides a 
sharp contrast between our approach and others that aim at weaker forms of 
independence between the paths, for example by seeking paths that are merely 
edge- or node-disjoint (UlIMlllZlIIlllllIllSllBllilJISS]. This is so because, as 
remarked elsewhere (e.g., [37]), independence by edge-disjointness encompasses 
independence by node-disjointness, which in turn encompasses independence by 
noninterference. Of course, the highest an independence relation's level in this 
hierarchy the easiest it is to implement it as the multiple paths are discovered 
(not coincidentally, the simple exchange of tokens between nodes suffices to 
produce edge- or node-disjoint path sets [29]). 

We proceed in the following manner. First we state the problem in graph- 
theoretic terms and give our solution in Section[2] Then we move, in Section|3l to 
a presentation of the methodology we followed in conducting our computational 
experiments. Our results are given in Section U] and involve comparisons with 
some prominent routing algorithms, viz. AODV [SF, AOMDV [31], OLSR [53], 
and MP-OLSR j54j. We used these algorithms both as stand-alone methods 
and as bases to our own heuristic. Our results include throughput and fairness 
[24] comparisons based both on NS2.34 [35] simulations and on the SERA link 
scheduling algorithm [48] . We continue in Section [5] with further discussions 
and conclude in Section |6j 

2 Problem formulation and heuristic 

For i and j any two distinct nodes of the WMN, we begin by assuming that 
some routing protocol already established a set Vij of paths directed from i 
to j. Another key element is that, since we seek to establish independence 
by noninterference, an interference model, along with its assumptions, must be 
selected. Our choice is the protocol-based interference model, together with the 
assumption that a node's communication and interference radii are the same. 



Should a different model be selected or the two radii be significantly different, 
the only effect would be for the graph construction process outlined below, once 
adapted accordingly, to produce a different graph (in particular, the interference 
radius could be chosen appropriately in order for the protocol-based interference 
model to mimic the physical interference model [42)'). 

Under the assumptions of the protocol-based interference model, the com- 
munication/interference radius is fixed at some value i?, which we take to be 
the same for all nodes. It follows that two nodes are neighbors of each other 
in the WMN if and only if the Euclidean distance between them is no greater 
than R. Moreover, since every link may transmit in both directions for error 
control, it also follows that two links can interfere with each other if either one's 
transmitter or receiver is a neighbor of the other's transmitter or receiver [8]. As 
will become apparent shortly, this has important implications when modeling 
interference, since links that share no nodes can still interfere with one another. 

In keeping with the locality principle outlined in Section [1] we work on the 
premise that a node fc's knowledge is limited to the set Nk of its neighbors and 
the set Nextfe(i, j) C Nk of neighbors to which it may forward packets sent by 
node i to node j (assuming k ^ j). Set Nextfc(i, j), obviously, depends on the 
j-bound paths that leave i and go through node k. The problem we study is 
that of eliminating from Vij the fewest possible paths (in a weighted sense, to 
be discussed later) so that the remaining path set, henceforth denoted by Vfj, 
contains no mutually interfering paths except at i or j. However, owing once 
again to the issue of locality, we forgo both optimality and feasibility a priori 
and settle for a heuristic instead. That is, for the sake of locality we admit 
the possibility that, in the end, neither will the selected paths be collectively 
optimal nor will the absence of interference among them be guaranteed (except 
at those paths' second hops, which will be non-interfering relative to one another 
necessarily) . 

Before proceeding, we tackle two special cases. The first one is that in which 
Vij contains the single- link path that connects i directly to j. In this case, we 
let V^ be the singleton that contains that path, since no other arrangement can 
possibly do better. The second special case is that in which \Vij\ = 1, provided 
the single path contained in Vij has at least two links, and is motivated by the 
situations in which Vij originates from single-path routing. In this case, we 
enlarge Vij before feeding it to our path-selection heuristic. Letting node k be 
such that Nexti(z,j) = {fc}, we do this enlargement of Vij by including in it 
as many paths from Vkj as possible, each suitably prefixed by the new origin 
i, provided none of the new paths is weightier than the one initially in Vij. If 
enlargement turns out to be impossible, then the problem becomes moot for Vij 
and we let Vf-=V,j. 

We are then in position to introduce our refinement algorithm for multi-path 
routing, henceforth referred to by the acronym MRA (for multi-path refinement 
algorithm) . The goal of MRA is to create an undirected graph Gij corresponding 
to the path set Vij and to extract from it the information necessary to determine 
V^. This graph's node set, henceforth denoted by V, has one node for each of 
the paths in Vij. Its edge set, denoted by E, is constructed in such a way as to 



represent every interference possibility that can be inferred solely from the sets 
Nk and Nextfc('i, j) for every k € Ni that participates in at least one of the paths 
in Vij. Once graph Gij is built, finding a maximum weighted independent set 
in it (i.e., a subset of V containing no two nodes joined by an edge in E and 
being as weighty as possible) provides the best possible approximate decision 
on which of the paths in Vij should constitute V^, namely those corresponding 
to the nodes in the maximum weighted independent set that was found. 

MRA proceeds as given next, following the introduction of some auxiliary 
nomenclature. Given two neighboring WMN nodes, say k and k' , we are inter- 
ested in the following two possibilities for the pair k,k' . A type- A pair is part 
of at least one path in Pij in such a way that k' £ Nextfc('i, j) with fc e iV^ or 
k G Nextfe'(i,j) with k' e Ni. In a type-B pair, both k and k' belong to at 
least one path in Vij as well, but not the same path. Moreover, at least one of 
them is a member of Ni while the other, if not in Ni as well, is in Nexti{i,j) for 
some other I £ N. Note that type- A and -B pairs constitute all the structural 
information that node i can gather by strictly local communication from its 
neighborhood Ni in the WMN. In the case of a type-A pair, we let Paths(fc, k') 
be the set of i-to-j paths to which the pair belongs. 

(1) Let V have one node for each path in Vij. Add one further node to V for 
each type-B pair of neighboring WMN nodes. We refer to these additional 
nodes as temporary nodes. 

(2) Construct E as follows: 

i. Let k, k' and I, V be two pairs of neighboring WMN nodes, each pair 
being of type A or B. If it holds that k = I, k = I', k' = I, or k' = I', 
then add an edge to E between each node corresponding to a path in 
Paths(fc, k') (if the pair is of type A) or the corresponding temporary 
node (if the pair is of type B) to each node corresponding to a path in 
Paths(l, I') (if the pair is of type A) or the corresponding temporary 
node (if the pair is of type B). 

ii. Connect any two nodes in V by an edge if, after the previous step, 
the distance between them is 2. 

iii. Remove all temporary nodes from V and all edges that touch them 
from E. 

(3) Find a maximum weighted independent set of Gij and output Vj^ accord- 
ingly 

In these steps, graph Gij starts out as a graph with \Vij\ nodes and gets 
enlarged by the addition of temporary nodes that represent some of the pos- 
sibilities of ofF-path interference as nodes in Ni engage in transmitting packets 
(Step (1)). Then it receives edges to account for the assumptions of the protocol- 
based interference model (Steps (2).i and (2).ii) and is after that stripped of all 
temporary nodes to end up with \Vij\ nodes once again (Step (2). iii). Step (2).ii, 
in particular, accounts for interference in the WMN when a link's transmitter or 



receiver does not coincide with (but is a neighbor of) another link's transmitter 
or receiver. The last MRA step, Step (3), is the determination of a maximum 
weighted independent set of Gij. We use node weights such that, for the node 
corresponding to path p g Vij, the weight is 1/Cp, where Cp is the path's hop 
count. In other words, shorter paths tend to be favored over longer ones as V^ 
is extracted from Vij. Clearly, though, any other desired criterion can be used 
as well. An illustration of how MRA works is given in Fig. [1] 

3 Methods 

We evaluated the performance of MRA through extensive experimentation with 
the following routing algorithms: AODV [38], AOMDV [3T], OLSR [23], and 
MP-OLSR [54. For the purpose of conciseness, we henceforth refer to the 
combination of each of these algorithms with MRA as R- AODV, R- AOMDV, R- 
OLSR, and R-MP-OLSR, respectively. We remark that, since both AODV and 
OLSR are single-path routing algorithms, handling the paths they generate falls 
into one of the special cases discussed in Section [2] Our experiments were run 
in the network simulator NS2.34 (NS2 henceforth) [35] and in a simulator that 
employs the SERA link scheduling algorithm [35] , briefly described later in this 
section. We used two different configurations of routing-algorithm parameters, 
and likewise two different configurations of NS2 parameters (one for the path- 
discovery process and another for performance evaluation) . These configurations 
were selected during initial tuning experiments and will be presented shortly. 

3.1 Topology generation 

We generated four types of network according to the maximum number of neigh- 
bors, A, a node may have. Each network was generated by placing n nodes inside 
a square of side 1500. The first of these nodes was positioned at the center of 
the square and the remaining nodes were placed randomly as a function of the 
communication/interference radius R introduced at the beginning of Section |2| 
Their placement was subject to the constraints that each node would have at 
least one neighbor, that no node would be closer to any other than 25 units 
of Euclidean distance, and that no node would have more than A neighbors. 
No more than 1000 attempts at positioning nodes were allowed; if this limit 
was reached then the growing network was discarded and the generation of a 
new one was started. The value of R was determined so that the expected den- 
sity of nodes inside a radius- i? circle would be proportional to A/i?^ (assuming 
some uniformly random form of placement), and be moreover about the same 
density as that of the whole network. It follows that A/i?^ oc n. We chose 
the proportionality constant to yield R = 200 for ri = 80 and A = 4, whence 
R = 200^20A/n. Of all the networks generated, there are 100 networks for 
each combination of n e {60, 80, 100, 120} and A G {4, 8, 16, 32}, thus totaling 
1600 networks. 
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Figure 1: Construction of graph Gij for the path set Vij ~ {a, 6,c, d, e}. The 
paths in Vij are shown in panel (a), where each directed edge leads to a node 
in some Next(i,j) set (for example, Nextfcj(i,j) = {j, ^5}) and each dashed 
edge joins neighboring nodes that share none of the i-to-j paths. We see in 
this panel that Paths(A;2, fee) = {b, c}, and also that there are five type-A node 
pairs (fci,j; fc2,fc6; fc2,fc3; ^1,^5; and k4,ke) as well as one single type-B node 
pair {k^jk^, labeled x). Panel (b) shows graph Gij as it stands after Steps (1) 
and (2).i, respectively for creating its node set V and initializing its edge set E 
as a function of node coincidence among all type-A or -B pairs. The distance- 
2 closure determined in Step (2).ii is shown in panel (c), with the additional 
edges represented as dashed lines. Panel (d), finally, shows Gij as it stands at 
the end, after Step (2).iii. Notice, in particular, the fundamental role played 
by node x in identifying the interference between paths c and d under the 
protocol-based interference model (it is through x, as the distance-2 closure is 
determined, that c and d become connected). If equal weights were used for all 
nodes, then clearly any of {a, 6}, {a, c}, {a, e}, {b, d}, or {d, e} would qualify as 
a maximum weighted independent set of Gij in Step (3). However, using 1/Cp 
as the weight of the node corresponding to path p leads to {a, b} and {a, e} as 
the only possibilities, each of total weight 1/2 + 1/3 = 5/6, since Cq = 2 and 
Cb = Ce = 3. In the end, then, we have either Vjj — {a, b} or V^ = {a, e}, each 
set containing paths whose second hops do not interfere with one another, each 
set as large as possible in the weighted sense provided by the use of 1/Cp. 



3.2 Path discovery 

For each of the 1600 networks we randomly generated 100 sets of node pairs to 
function as origin-destination pairs (instances of the i,j pair we have been using 
throughout). Each set comprises n pairs and no node was ahowed to appear 
more than once in any set as an origin. For each of the node-pair sets and 
each of the four routing algorithms (AODV, AOMDV, OLSR, and MP-OLSR) 
we obtained n path sets (instances of Vij). Likewise, for each of the node- pair 
sets and each of the four refined routing algorithms (R-AODV, R- AOMDV, 
R-OLSR, and R-MP-OLSR) we obtained another n path sets (instances of V^). 

For each routing algorithm, the discovery of each Vij instance (i.e., for a 
single origin and a single destination) proceeded as follows. After loading the 
network topology onto NS2 we conducted a 15-second simulation with one flow 
agent for the single origin i and the single destination j, using a CBR of one 
1000-byte packet per second. The remaining pairs were handled likewise after 
resetting the simulator. We remark that this one-pair-at-a-time strategy, as 
opposed to generating paths for all n pairs in the same set concomitantly, was 
meant to minimize packet loss due to path overload and also to avoid the possible 
interference of a previously discovered path with the discovery of a new one. Of 
course, this argument is only valid for on-demand routing algorithms (AODV 
and its variants depend on network load, while the OLSR variants always find 
identical paths for the same network topology), but we proceeded in this way 
in all cases. As a consequence, our experiments are entirely reproducible. 

We set NS2 to its default configuration, but employed the DRAND MAC 
protocol [40] to avoid collisions in the path-discovery process. To adjust the 
radius R we also set the parameter RXThresh- (RXT) to the appropriate value 
given by the program threshold. cc (cf. the NS2 manual). We used the im- 
plementations of the AODV, OLSR, and AOMDV routing agents available in 
version 2.34 of NS2 and the MP-OLSR routing agent available at [32]. For 
AOMDV, we made the small modifications proposed by [5B] to discover only 
node-disjoint paths with at most K paths for each origin-destination pair of 
nodes. We adopted these modifications because node-disjoint paths are clearly 
more interference-free than otherwise. We chose K = 5 because it achieved the 
best throughput values for 2 < K < 7. The same modifications were effected on 
MP-OLSR (as proposed by [57]). Out of the same range for K, and for the same 
reason as above, we used ii' = 3 for A e {4, 8} and ii' = 5 for A £ {16, 32}. 

3.3 Performance evaluation 

Once we fix a value for n and a value for A, there are 10^ path sets on which 
to evaluate the performance of MRA. Each of these sets is relative to n origin- 
destination pairs. In our experiments, we randomly grouped these pairs into n 
sets, each containing a different number of origin-destination pairs (i.e., one set 
containing a single pair, another containing two pairs, and so on), and simulated 
the network's behavior in transporting predominantly heavy traffic from the 
origins to the destinations. We use OD to denote the set of pairs in question. 



Table 1: Parameters used for tuning. 



NS2 parameter 


Interval 


Increment 


CBR interval 0.001s . . 0.005s 0.0001 
CST 0.1RXT..2RXT O.IRXT 

Data/basic rate 1/1, 2/1, 11/2 Mbps 

Table 2: Parameters used for performance evaluation. 


NS2 parameter 




A 






4 


8 


16 32 



CBRl interval 0.0025s 0.0027s 0.0029s 0.0031s 

CBR2 interval 0.0045s 

CBR packet size 1000 bytes 

CST 0.6RXT 0.7RXT 0.8RXT 0.9RXT 
Data/basic rate 11/2 Mbps 

hence 1 < \0D\ < n, and let V = [J.jeOD "^y ^^d P^ ^ UijeOD ^5- For the 
sake of normalization, all our performance results are presented against the pair 
density = \OD\/n e (0,1]. 

Each experiment began by loading the corresponding network onto NS2, 
with MAC set to 802.11 and the routing agent to NOAH [34]. NOAH works 
only with fixed paths that have to be configured manually and therefore does not 
send routing-related packets, thus providing the ideal setting for a performance 
evaluation free of any interference from control packets. Next we started a CBR 
traffic flow from each origin to each destination and measured the number of 
successfully delivered packets during the last 120 of the 135 seconds of simulation 
(following, therefore, a warm-up period of 15 seconds). 

During the initial, tuning experiments we varied three of the NS2 parameters 
widely. These were the carrier sense threshold (CST), aiming to increase the 
spatial reuse and consequently the throughput [IS] , and the CBR parameter as 
well as the network transmission rate, aiming to obtain as much throughput and 
fairness as possible (cf. Section 14. 2p . Table [T] presents the parameter ranges of 
our tuning experiments, during which the highest CBR rate that afforded some 
gain in throughput was identified for each value of A, and also lowest rate that 
afforded some gain in fairness]^ Our choices for use thereafter while conducting 
performance evaluation are shown in Table [21 where CBRl and CBR2 refer to 
such highest and lowest rate, respectively. 

All our NS2 experiments were carried out under 802.11, which is a CSMA 
protocol. As an alternative setting that might provide some insight into the 
performance of MRA under some TDMA scheme (an approach fundamentally 
distinct from CSMA [IH]), we selected the SERA link scheduling algorithm [15] . 



^We observed in NS2 simulations that, within Umits, increasing the CBR rate^ parameter 
tends to lead to an increase in throughput while decreasing it tends to lead to an increase in 
fairness. 



SERA seeks to schedule the Unks of a set of paths while striving to maximize 
throughput on those paths. It is therefore quite well suited to the task at 
hand. The throughput that our SERA simulator provides is given in terms of 
time slots, so in order to achieve a meaningful basis for comparing CSMA- and 
TDMA-based results a translation is needed of such throughput figures into 
those provided by NS2 in the CSMA experiments. We did this by resorting 
to a very simple NS2 simulation to determine the duration of a time slot. In 
this experiment, a node sends packets to another and the time for successful 
deliveries is recorded. Since this time is conceptually the same that in SERA 
is taken to be a time slot, the translation from one setting to the other can be 
accomplished easily. Using the parameter values shown in Table [21 we found a 
time-slot duration of 0.002 seconds. This is the duration we use, together with 
the ND-BF numbering scheme for SERA and its B parameter set to 2 (cf. [48]). 

4 Computational results 

We divide our results into two categories. First we present a statistical analysis 
of the networks generated and the path sets obtained by the original routing 
algorithms and by their refinements through the use of MRA. Then we present 
the ratios of the refined algorithms' throughputs to those of their corresponding 
originals (absolute values are given in Figs. IA.1HA.3]) and also fairness figures. 
For the purpose of conciseness we report only on the n = 120 results, since they 
are qualitatively similar to those related to the other three values of n we used. 
In Section [5l though, we do discuss some of the quantitative differences that 
were observed. 

4.1 Properties of the networks generated 

The 1600 networks we generated are that same that were used in [3^. We re- 
fer the reader to Section 8.1 of that publication for a variety of the networks' 
statistical properties, such as the occurrence of topologies structured in some 
particular way and some of their structure-related distributions. Here we con- 
centrate on presenting those properties that pertain to routing both before and 
after refinement, since they are the ones we have found useful in helping explain 
the throughput results we present later. 

The average path multiplicity (number of paths) per origin-destination pair 
in the path sets V and V^ is given in Table [3] for every combination of A and 
\OD\. Note that MP-OLSR is absent from the table in spite of being a multi- 
path algorithm, the reason being that in this case the K parameter does not 
work as an upper bound (as it does for AOMDV), but rather as the fixed number 
of paths to be found. We observe in the table that the average path multiplicity 
increases monotonically with A for fixed lOUj, which is expected from the well- 
known fact that the number of possible paths between the same two nodes 
grows with A in arbitrary graphs [25]. On the other hand, increasing \0D\ for 
fixed A causes very little variation, probably owing to the method we used to 
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Table 3: Average path multiplicity per origin-destination pair for n = 120. Data 
are averages over 10* instances of V or V^ for each combination of A and \0D\. 



A 


OD\ 


R-AODV 


AOMDV 


R-AOMDV 


R-OLSR 


R-MP-OLSR 


4 


1..10 


1.5 


3.3 


1.7 


1.4 


1.5 




11.. 90 


1.5 


3.4 


1.8 


1.5 


1.6 




91 . . 120 


1.6 


3.4 


1.8 


1.5 


1.6 


8 


1..10 


1.6 


3.6 


1.9 


1.6 


1.7 




11.. 20 


1.7 


3.6 


2.1 


1.6 


1.7 




21 . . 70 


1.7 


3.7 


2.1 


1.6 


1.7 




71 . . 80 


1.7 


3.7 


2.1 


1.6 


1.8 




81 . . 120 


1.7 


3.7 


2.2 


1.6 


1.8 


16 


1..10 


2 


3.9 


2.9 


2 


2 




11.. 50 


2.1 


3.9 


2.9 


2 


2 




51.. 90 


2.1 


4 


3 


2.1 


2 




91 . . 120 


2.1 


4 


3 


2.1 


2.1 


32 


1..10 


2.2 


4.1 


3.1 


2.2 


2.2 




11.. 60 


2.2 


4.1 


3.1 


2.2 


2.3 




61 . . 70 


2.3 


4.2 


3.2 


2.2 


2.4 




71 . . 80 


2.3 


4.2 


3.2 


2.2 


2.5 




81.. 90 


2.3 


4.2 


3.2 


2.2 


2.6 




91 . . 120 


2.3 


4.2 


3.2 


2.3 


2.7 



discover the paths in the first place (i.e., by handling each origin-destination 
pair independently of the others). The overall pruning effect of MRA is also 
clearly manifest in the table, since in all cases the average for an algorithm's 
refined version is less than that of the original algorithm (i.e., on average we 
have \V'^\ < |7^| for all routing algorithms). 

The path-size distributions for V and V^ with \OD\ = n (that is, for the 
original algorithms and their refinements, using the full path sets, as generated) 
are given in Figs. [2] and [S] respectively, for every value of A. Note first that, as 
expected by virtue of the well-known dependency of path sizes on the average 
degree of nodes [17], increasing A leads the path-size distribution for a given 
algorithm to peak at ever smaller values. Another expected result is that the 
OLSR variants all allow for longer paths than those of AODV. The reason for 
this is that the multi-point relay (MPR) concept at the heart of OLSR tends 
to cause longer paths to be produced than the shortest-path algorithm (SPA) 
used by AODV. It is also noteworthy that the application of MRA to yield the 
refined algorithms has no impact on the distributions other than causing some 
of the greatest path sizes to occur more frequently. 

4.2 Results 

Our computational results relating to throughput are summarized in Fig. 2] for 
CBRl, Fig. [5] for CBR2, and Fig. [6] for SERA. Each figure is organized as 
a set of four panels, each for one of the four values of A. All three figures 
show the behavior of the ratio, here denoted by a, of each refined algorithm's 
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Figure 2: Distribution of path sizes for the original algorithms with n = 120, 
for A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). For each value of 
A the distribution refers to 100 networks and 100 path sets per network, each 
corresponding to n origin-destination pairs. 
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Figure 3: Distribution of path sizes for the refined algorithms with n = 120, 
for A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). For each value of 
A the distribution refers to 100 networks and 100 path sets per network, each 
corresponding to n origin-destination pairs. 
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throughput to that of its original version. All plots are given against the ratio 
9 introduced earlier, which indicates what fraction of the n origin-destination 
pairs corresponding to each path set is being taken into account. In all panels 
a highlighted horizontal line is used to mark the a = 1 threshold, above which 
a refined algorithm can be said to have overtaken its original version. 

Except for a few cases of < 0.1 (in which only a few origin-destination pairs 
coexist in the network and therefore refinement to achieve path-independence 
by noninterference is probably pointless to begin with), it follows from Figs. 3]- 
[S] that MRA was effective to some extent in all cases. Fixing the simulation 
scenario (CBRl, CBR2, or SERA) reveals that the behavior of a depends very 
little on the value of A or 6* (provided 9 is sufficiently large). Overall, the 
values of a seem best for the coupling of OLSR with SERA, followed by CBRl 
coupled with either OLSR or MP-OLSR, and lastly for CBR2 without any 
marked preference for any routing method. Recall that the SERA scenario 
is TDMA-based, while both CBRl and CBR2 are CSMA-based, with CBRl 
operating at the higher rates. 

Another perspective from which it is worth examining performance is that 
of fairness in the distribution of traffic through the paths. In other words, 
given an origin-destination pair and the multiple paths leading from the origin 
to the destination, we look at how traffic gets distributed through the various 
paths. One way of quantifying this is by means of the fairness index pi] . 
Given a set of paths Qo the corresponding fairness index can be defined as 
(SpGQ ^p)^/|2l SpsQ -^p' '^tisre Xp is the number of packets delivered to j 
through path p during the experiment. The fairness index ranges from 1/|Q| to 
1, indicating when equal to 1 that traffic is evenly distributed among the paths. 

We give results on the fairness index in Figs. [7HS1 respectively for CBRl, 
CBR2, and SERA. Note initially that, somewhat unexpectedly (owing to the 
algorithms' markedly different strategies), AODV and OLSR are statistically 
indistinguishable from each other as far as fairness is concerned. The same holds 
for their refined versions, respectively R-AODV and R-OLSR. Note also that 
the fairness index, in all cases, tends to decrease as 9 is increased. This means 
that, as might be expected, the presence of denser end-to-end traffic tends to 
disrupt the balance between paths more easily. AODV and OLSR have the best 
figures overall, better even than their refined versions. So, unlike throughput, 
fairness does not seem to improve as we extend a single-path algorithm's set of 
paths and then apply MRA for refinement. All multi-path strategies, on the 
other hand, can be seen to benefit from the use of MRA. 

Another interesting trend that can be observed in Figs.[7HS]is that, though 
always decaying with 9 as we noted, in general the fairness index is best for 
SERA, followed by CBR2, then by CBRl. While we believe the position of 
SERA in this rank to be closely related to its TDMA-based nature, it is curious 
to observe that the relative positions of CBRl and CBR2 are exchanged with 
respect to what we observed for throughput. It seems, then, that in selecting 



■^We use Q as a place holder for either V or V^, depending on whether the routing algorithm 
in question is one of the originals or one of the refinements through MRA. 
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Figure 4: Throughput ratio for CBRl with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 0. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure 5: Throughput ratio for CBR2 with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 0. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure 6: Throughput ratio for SERA with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 0. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure 7: Fairness index for CBRl with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 0. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure 8: Fairness index for CBR2 with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 9. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure 9: Fairness index for SERA with n = 120, for A = 4 (a), A = 8 (b), 
A = 16 (c), and A = 32 (d). Data are averages over the 10* path sets that 
correspond to each value of A for each value of 0. Confidence intervals are less 
than 1% of the mean at the 95% level, so error bars are omitted. 
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between the rates associated with CBRl and CBR2 one is automaticahy forced 
to favor throughput over fairness or conversely. 

Our definition of the fairness index given above is only one of the possibilities 
in the context of multi-path routing. Another alternative is to coalesce all pack- 
ets delivered from one origin, say i, to one destination, say j, into a single num- 
ber Xij and then compute the fairness index as (J^ijeOD ^ij)'^/\0^\ J2ijeOD ^Ij- 
Proceeding in this way would shift the focus of the fairness index from paths 
to origin-destination pairs. We give no detailed results on this alternative, but 
do provide an example with SERA in Fig. IA.4I There is clearly great similarity 
with the results in Fig. [3J but the numbers tend to be higher. 



5 Discussion 

As we remarked earlier, we have provided no results for n G {60, 80, 100} be- 
cause, essentially, they are indistinguishable from those of n = 120 in qualitative 
terms. We do remark, however, that a few noteworthy quantitative differences 
were observed. For example, in the case of SERA a higher throughput ratio a 
was sometimes observed for reduced n but the same value of the pair density 
9. This stems not only from the fact that the average path size decreases as 
n decreases, but more generally from the fact that the number of links for the 
same 9 is lower in the smaller networks, thus fewer links interfere with one an- 
other and fewer links have to be scheduled. Such improvement in the value of 
a, therefore, seems to depend on the network's path-size distribution. 

As we noted briefly in Section H^l often OLSR and MP-OLSR turned out 
to be the routing methods most prone to benefit from the refinement provided 
by MR A. One clue as to why this is so may be already present in Table [31 
where the refined versions of these two methods have some of the lowest path 
multiplicities overall, thence a tendency to incur less interference. However, 
this holds also for the refined version of AODV, so there has to be some other 
distinguishing aspect. While our data are not sufficient to provide a definitive 
answer, we believe the two methods' superiority to be owed to a combination of 
MRA (which attempts to reduce path multiplicity to lower interference) with 
the MPR concept that is intrinsic to the OLSR variants (which in turn attempts 
to provide an initial set of paths as spatially distributed as possible). As for 
AODV, the SPA at its core probably produces paths that are less spatially 
separated [551 [TU] . 

In a related vein, the results of Section W^ also point at SERA as providing 
superior throughput results vis-a-vis those of CBRl and CBR2, and similarly 
for CBRl with respect to CBR2. Because such results refer to throughput 
gains after refinement by MRA, more detailed data are needed for a direct 
comparison of throughputs. These are shown in Figs. IA.5HA.7l which clearly 
confirm the rank. Once again the reason for this is not totally clear, but it 
may be a manifestation of SERA's properties rather than of the superiority of 
its underlying TDMA scheme over the CSMA protocols of CBRl and CBR2. 
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After all, a considerable amount of research has been directed toward the TDMA 
versus CSMA question, without however reaching an agreement [161 HOl HSl E] ■ 
Another curiosity related to this issue of how SERA compares to CBRl and 
CBR2 has to do with the total absence of traffic on some paths. While SERA 
precludes this from occurring as a matter of design principle, unused paths 
do occur in the other two cases. During the corresponding NS2 simulations 
this persisted even if simulation times were extended or the synchronization of 
multiple CBRs was reordered or entirely removedjfl As it turns out, it seems 
that certain paths remained unused so that throughput could be increased on 
other paths. Data on the most critical cases, viz. those in which no packets at 
all were delivered for an origin-destination pair, are shown in Figs . IXfSl and [XTOl 
Examining these data confirms our expectation that the single-path algorithms 
should be less prone to the occurrence of such extreme cases. It also reveals 
that both R-AOMDV and R-MP-OLSR had fewer such occurrences than their 
corresponding originals (i.e., MRA seems to have attenuated the problem). 

6 Concluding remarks 

MRA is a heuristic for the refinement of routing paths in WMNs. It was devel- 
oped with multi-path routing algorithms in mind but works also on single-path 
algorithms (through a manipulation of its input to obtain multiple paths from 
the overall set of single paths). MRA is fully local, in the sense that it depends 
only on information that is readily available to each node and its immediate 
neighborhood in the network. Being local means that the extra control traffic it 
may entail is negligible, especially if we consider that it need happen only once 
for a fixed set of paths. It involves the solution of an NP-hard problem, that 
of finding a maximum weighted independent set in a graph, but the inputs in- 
volved are typically very small, leading to negligible running times if compared 
to arbitrary graphs |36) . 

Our computational experiments with MRA on a TDMA setting (SERA) and 
two CSMA settings (CBRl and CBR2) revealed improvements in throughput 
of up to 30% for both the AODV and OLSR routing methods and their multi- 
path variants. The latter were also improved by MRA in terms of fairness. 
Finding out just how close these improvements come to the very best that can 
be achieved is an open problem that involves solving the path-pruning problem 
globally. This problem is NP-hard as the one solved by MRA, but its input, 
relating to the WMN in a global scale, is much more sizable. 

Another interesting aspect for further investigation is the effect of MRA on 
multi-radio networks. In such schemes the network's path density can be dras- 
tically reduced for a given frequency, thus potentially benefiting MRA. Con- 
versely, it may also be possible to use MRA to achieve the desirable goal of 
minimizing the number of radios [5J 139) . Further research is also needed on the 
trade-off that clearly exists between obtaining spatially separated paths or rel- 



^As these appear to be commonly occurring problems of NS2 agents. 
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atively short ones. While the former is good for non-interference, disregarding 
the latter may lead to throughput loss on the relatively longer paths. 
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Figure A.l: Throughput (packets/time slot) for CBRl with n = 120, for A = 4 
(a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 10* path 
sets that correspond to each value of A for each value of 0. Confidence intervals 
are less than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure A. 2: Throughput (packets/time slot) for CBR2 with n = 120, for A = 4 
(a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 10* path 
sets that correspond to each value of A for each value of 9. Confidence intervals 
are less than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure A.3: Throughput (packets/time slot) for SERA with n = 120, for A = 4 
(a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 10^ 
path sets that correspond to each value of A for each value of Q. Error bars are 
based on confidence intervals at the 95% level. 
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Figure A. 4: Alternative fairness index for SERA with n = 120, for A = 4 (a), 
A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 10^ path 
sets that correspond to each value of A for each value of 9. Confidence intervals 
are less than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure A. 5: Ratio of SERA's throughput to that of CBRl with n = 120, for 
A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 
10^ path sets that correspond to each value of A for each value of 6. Confidence 
intervals are less than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure A.6: Ratio of SERA's throughput to that of CBR2 with n = 120, for 
A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 
10^ path sets that correspond to each value of A for each value of 6. Confidence 
intervals are less than 1% of the mean at the 95% level, so error bars are omitted. 
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Figure A.7: Ratio of CBRl's throughput to that of CBR2 with n = 120, for 
A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 (d). Data are averages over the 
10^ path sets that correspond to each value of A for each value of Q. Confidence 
intervals are less than 1% of the mean at the 95% level, so error bars are omitted. 



35 



(a) 
0.50 

0.40 
§ 0.30 
& 0.20 

0.10 

(b) 
0.50 

0.40 
§ 0.30 
ffi 0.20 

0.10 

(c) 
0.50 

0.40 

§ 0.30 

[i; 0.20 

0.10 

(d) 
0.50 

0.40 

§ 0.30 

I 

i: 0.20 

0.10 
0.00 



AODV ^»»o»>a 

R-AODV ^^ 

AOMDV ^^ 

R-AOMDV ^ 

OLSR N\\\\\l 

R-OLSR 

MP-OLSR vm////m 

R-MP-OLSR ^^ 



^a 




AODV txwxwj 

R-AODV ^^ 

AOMDV ^^ 

R-AOMDV ^ 

OLSR N\\\\\1 

R-OLSR ^^ 

MP-OLSR tM-////////a 

R-MP-OLSR ^^ 



AODV tooxf/i 

R-AODV ^^ 

AOMDV 

R-AOMDV 

OLSR NWWM 
R-OLSR 
MP-OLSR ymm/A 
R-MP-OLSR ^^ 



AODV txx^KXi 
R-AODV ^^ 
AOMDV Bww^ 

R-AOMDV I 1 

OLSR NWWM 

R-OLSR ^^ 

MP-OLSR '«"'«""'«'. 

R-MP-OLSR ^^ 



^ 1 Ji I ^ -I ti 



j_i 






15 20 25 30 40 50 60 70 80 90 100 110 120 
Number of origin-destination pairs 



Figure A. 8: Distribution of the number of origin-destination pairs without any 
traffic in CBRl with n = 120, for A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 
(d). Data are averages over the 10"* path sets that correspond to each value of 
A. 
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Figure A. 9: Distribution of the number of origin-destination pairs without any 
traffic in CBR2 with n = 120, for A = 4 (a), A = 8 (b), A = 16 (c), and A = 32 
(d). Data are averages over the 10"* path sets that correspond to each value of 
A. 
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